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O I This presentation reports on recent developments concerning basic aspects of low-energy 

^ I QCD as they relate to the understanding of the nucleon mass and the nuclear many-body 

y—t • problem. 
> '. 

in . 

l> . 1. INTRODUCTION: SOME BASICS OF LOW-ENERGY QCD 
O ; 

^ . Nuclear physics has a well developed phenomenology in terms of nucleons interacting 
, through potentials constrained by nucleon-nucleon scattering data. At the same time, 
^ I nucleons and nuclei are aggregates of quarks and gluons residing in the hadronic, low- 
I temperature phase of QCD. In recent years, the understanding of the nucleon itself, 
^ ■ in particular of its low-energy, long-wavelength structure, has been progressing rapidly. 
^ . Likewise, steps have been taken to explore the possible connection of the nuclear many- 
P3 ! body problem with low-energy QCD, its symmetries and symmetry breaking patterns. It 
j> I is therefore timely and appropriate to discuss these two related themes from a unifying 
I perspective. 

I 1.1. Spontaneous chiral symmetry breaking and the hadronic mass gap 

At low energies confinement implies that the QCD eigenstates are not elementary quarks 
and gluons but hadrons, colour singlet composites of quarks and gluons. Confinement also 
implies that the chiral SU{2)l x SU{2)r symmetry of two-flavour QCD with (almost) 
massless u- and d-quarks is spontaneously broken. As a consequence, the QCD vacuum 
is non-trivial: it hosts a strong condensate of (scalar) quark-antiquark pairs, the chiral 
condensate (gg) ~ — 1.5fm^^. Goldstone's theorem states that there exist low frequency 
collective excitations of the condensed QCD ground state: pseudoscalar Goldstone bosons, 
identified with the pions. Their mass vanishes in the exact chiral limit, i.e. for massless 
u- and (i-quarks. The finite but small physical pion mass, uIt^ ~ 0.14 GeV, reflects the 
non-zero quark masses rrtu^ ~ 5 MeV. Goldstone's theorem also states that these pions 
interact weakly at low energy and momentum. In fact, their interaction strictly vanishes 
in the exact chiral limit and at zero energy and momentum. 

A charged pion decays through weak interaction into a muon and an antineutrino. 
The hadronic matrix element of the axial vector current connecting the pion with the 
QCD vacuum determines the pion decay constant, with the empirical value /jr = 
92.4 ± 0.3 MeV. This constant reflects, as an order parameter, the spontaneously broken 
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chiral symmetry of the QCD ground state just like the chiral (quark) condensate itself. 
The symmetry breaking pattern of low-energy QCD is manifest in a characteristic mass 
gap that governs the spectrum of the lightest hadrons. The scale associated with this 
mass gap is A ~ 47?/,^ ~ 1 GeV. 

Low-energy QCD is the physics of strong interactions at energies and momenta small 
compared to the 1 GeV scale of spontaneous chiral symmetry breaking. In this low- 
frequency, long- wavelength limit QCD is realized in the form of an effective field theory 
of weakly interacting Goldstone bosons coupled to heavier hadrons which act as almost 
static sources. Given this scenario, we address the following questions: 

• How does the nucleon develop its mass, starting from almost massless quarks, and 
what role does the pion cloud of the nucleon play in this context? 

• What is the impact of spontaneously broken chiral symmetry in QCD on the nuclear 
many-body problem? 

1.2. Chiral effective field theory 

The mass scale set by the gap A ~ 47r/^ offers a natural separation between "light" 
and "heavy" (or, correspondingly, "fast" and "slow") degrees of freedom. The basic 
idea of an effective field theory is to introduce the active light particles as collective 
degrees of freedom, while the heavy particles are treated as almost static. In low-energy 
QCD, the "light" particles are the pions whereas the nucleons figure as "heavy". The 
low-energy strong interaction dynamics is described by an effective Lagrangian which 
incorporates all relevant symmetries of QCD. The elementary quarks and gluons of the 
original Lagrangian are now replaced by the Goldstone bosons of spontaneously broken 
chiral symmetry. Heavy particles, such as nucleons, act as sources for the Goldstone 
bosons. Their leading-order couplings are determined by the Noether currents of the 
underlying symmetries. 

In the sector of QCD with the two lightest quark flavours, relevant to most of nuclear 
physics, the Goldstone bosons are represented by a 2 x 2 matrix field U (x) G SU{2) which 
collects the three isospin components iTaix) of the pion. A convenient representation is 
U{x) = exp[iTa7^a{x)/ f] where the pion decay constant / in the chiral limit provides 
a suitable normalisation. In essence, the QCD Lagrangian is replaced by an effective 
Lagrangian expressed in terms of U{x) and its derivatives: 



Goldstone bosons do not interact unless they have non-zero four-momentum, so the low- 
energy expansion of Ceff is an ordering in powers of (9^t/. 

The effective Lagrangian in the pure meson sector (with zero baryon number) has the 
following leading terms [ilj: 



The second term on the r.h.s. represents the (perturbative) explicit chiral symmetry 
breaking by the quark mass matrix, m = diag{mu,rnd)- The constant B is related to 
the chiral condensate as Bp = —{qq) = rn^p/imu + m^). The last equality follows 



CQCD^c,ff{u,du,d^u,...). 



(1) 



Ceff = C^= '-Tr[d^U^d^U] + '-BTr[m{U^ + U)] + 0[{dUY] . 



(2) 



3 



from identifying the pion mass term in the Lagrangian. With inclusion of higher order 
corrections the systematic expansion of the S matrix based on Chiral Perturbation 
Theory (ChPT), has become a powerful tool to deal with low-energy observables. One of 
its prime successes is the quantitative description of pion-pion scattering close to threshold 

m 

Consider next the sector with one unit of baryon number which includes the physics 
of the interacting pion-nucleon system. We restrict ourselves to the case of Nf = 2 
flavours. The previous pure meson Lagrangian is now replaced by Ceff = '^■k + '^■kN- The 
additional term involving the nucleon is expanded again in powers of derivatives (external 
pion momenta) and quark masses. The vrA^ effective Lagrangian to second order in the 
pion field has the form 

CnN = ^'^(^7/.^'' - M;v)^iv - ^tvTmTsT^tv- a'^TT (3) 



The Dirac spinor field = (p? of the nucleon is represented as an isospin-1/2 doublet 
of proton and neutron. The nucleon mass 

M^ = Mo + a^ , (4) 

has a large part Mq which exists already in the chiral limit of vanishing quark masses. 
The sigma term 

o-Tv = {N\mg{uu + dd)\N) (5) 

is the correction to Mq from the small quark mass term of the QCD Lagrangian, with 
niq = {rrtu + md)/2. Its empirical value is — 45 — 55 MeV [Oj. The vector and 
axial vector couplings of the nucleon to the Goldstone boson fields are dictated by chiral 
symmetry. They involve known structure constants: the pion decay constant and the 
nucleon axial vector coupling constant, qa = 1.270 ± 0.003, determined from neutron 
beta decay. Not shown in Eq.(3) is a series of additional terms of order [d^n)"^. These 
terms come with further constants that need to be fitted to experimental data and reflect, 
for instance, the very important role played by the A(1232) resonance in pion-nucleon 
scattering. 

Baryon ChPT, the calculational framework based on C^ff = '^n + '^nN, has been applied 
quite successfully to a variety of low-energy processes (such as pion-nucleon scattering, 
threshold pion photo- and electroproduction and Compton scattering on the nucleon) for 
which increasingly accurate experimental data have become available in the last decade. 
Reviews can be found in [ 4]. In recent years the nucleon-nucleon interaction at long and 
intermediate distance scales has also become an active area of ChPT [ . 

2. THE NUCLEON: ITS MASS AND SCALAR FIELD 

2.1. Nucleon mass and pion cloud 

Understanding the nucleon mass is clearly one of the most fundamental issues in nuclear 
and particle physics [ [7j. Progress is now being made towards a synthesis of lattice 
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QCD and chiral effective field theory, such that extrapolations of lattice results to actual 
observables are beginning to be feasible [|HIE]- Accurate computations of the nucleon mass 
on the lattice have become available [ ^1 El E] , but so far with u and d quark masses 
exceeding their commonly accepted small values by typically an order of magnitude. 
Methods based on chiral effective theory can then be used, within limits, to interpolate 
between lattice results and physical observables. 

The nucleon mass is determined by the expectation value (iV|0|^|iV) of the trace of the 
QCD energy-momentum tensor, 9^ = {P{g)/2g)G^^G^''' + muUU + mddd+ ... , where G^" 
is the gluonic field tensor, P{g) is the beta function of QCD, and qq with q = u,d, ... 
are the quark mass terms (omitting here the anomalous dimension of the mass operator 
for brevity). Neglecting small contributions from heavy quarks, the nucleon mass (4) 
taken in the SU(2)f chiral limit, ^ — >■ 0, is 

Mo = {N\fG,,G^^\N) . (6) 

This relation emphasises the gluonic origin of the bulk part of M^r, the part for which 
lattice QCD provides an approriate tool to generate the relevant gluon field configurations. 
At the same time, QCD sum rules connect Mq to the chiral condensate (qq). The leading- 
order result (loffe's formula [^|) is: 

Mo = -^{qq) + - , (7) 

where ~ 1 GeV is an auxiliary scale (the Borel mass) which separates "short" and 
"long" distance physics in the QCD sum rule analysis. While loffe's formula needs to be 
improved by including condensates of higher dimension, it nevertheless demonstrates the 
close connection between dynamical mass generation and spontaneous chiral symmetry 
breaking in QCD. For nuclear physics, Eq.(7) gives important hints [^5: the change 
of the condensate (qq) with increasing baryon density implies a significant reduction of 
the nucleon mass in the nuclear medium. This observation finds its correspondence in 
the strong scalar field as it appears in phenomenological relativistic mean field models [ 

caiinj. 

In chiral effective field theory, the quark mass dependence of M^v translates into a 
dependence on the pion mass, ~ m^, at leading order. The dressing of the nucleon 
with its pion cloud, at one-loop order, is illustrated in Fig.l. The systematic chiral 
expansion of the nucleon mass gives an expression of the form [ 9 : 

= Mo + cml + dmt - ^-^g\m. ^ " ^) + ' 

where the coefficients c and d multiplying even powers of the pion mass include low-energy 
constants constrained by pion-nucleon scattering. The coefficient d also involves a log m.j^ 
term. Note that the piece of order (non-analytic in the quark mass) is given model- 
independently in terms of the known weak decay constants gA and /,r (strictly speaking: 
by their values in the chiral limit). 

The interpolation shown in Fig. 2 determines the nucleon mass in the chiral limit, Mq, 
and the sigma term ctat = J2q=u,d^q{dMN / druq) . One finds Mq ~ 0.89 GeV and aN = 
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(47±3) MeV [Ej in this approach. (Treating the A isobar as an exphcit degree of freedom 
rather than absorbing it in low-energy constants, there is a tendency for shghtly larger 
values of (Jat). 
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Figure 1. Pion cloud contributions to 
the nucleon mass generated by the chi- 
ral effective Lagrangian (3). 



Figure 2. Best fit (solid curve) interpolating be- 
tween lattice results [ HUl CH 112] and the phys- 
ical nucleon mass, using NNLO chiral pertur- 
bation theory [ 9j. The dashed and dash-dotted 
curves show consecutive steps in the expansion [ 



2.2. Scalar form factor of the nucleon 

The prominent role played by the pion as a Goldstone boson of spontaneously broken 
chiral symmetry has its strong impact on the low-energy structure and dynamics of nu- 
cleons. When probing the individual nucleon with long-wavelength fields, a substantial 
part of the response comes from the pion cloud, the "soft" surface of the nucleon. While 
these features are well known and established for the electromagnetic form factors of the 
nucleon, its scalar-isoscalar meson cloud is less familiar and frequently obscured by the 
notion of an "effective sigma meson". On the other hand, the scalar field of the nucleon is 
at the origin of the intermediate range nucleon-nucleon force, the source of attraction that 
binds nuclei. Let us therefore have a closer look, guided by chiral effective field theory. 

Consider the nucleon form factor related to the scalar-isoscalar quark density, Gsiq"^) = 
{N{p')\uu + dd\N{p)), at squared momentum transfer = {p — p'Y- In fact, a better 
quantity to work with is the form factor (JNio^) = ^qGs{<f) associated with the scale 
invariant object mq{uu + dd). Assume that this form factor can be written as a subtracted 
dispersion relation: 

71 JAml t[t + Q^) 

where the sigma term aN introduced previously enters as a subtraction constant. We are 
interested in spacelike momentum transfers with = —q^ > 0. The dispersion integral 
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in Eq.(19) starts out at the two-pion threshold. It involves the spectral function rjsit) 
which includes all J'^ = O"*", / = excitations coupled to the nucleon: a continuum of even 
numbers of pions added to and interacting with the nucleon core. 




Figure 3. (a) Sketch of 
the scalar form factor of the 
nucleon; (b) Two-pion ex- 
change interaction between 
nucleons. 



Figure 4. Spectral function of the scalar-isoscalar 
nucleon form factor: chiral perturbation theory 
(ChPT at 2 loop NNLO [HT], solid and dashed) 
in comparison with empirically deduced spectrum 
[ ITH] (upper curve). The separate contribution of 
the nucleon Born term is also shown (dash-dotted). 



Chiral perturbation theory at next-to-next-to-leading order (NNLO) in two-loop ap- 
proximation has been applied in a recent calculation of the spectral function ris{t) [I17j. 
This calculation includes not only nucleon Born terms and leading vrvr interactions but 
also important effects of intermediate A isobar excitations in the two-pion dressing of the 
nucleon core (see Fig.Sa). The result (Fig.4) can be compared with the "empirical" scalar- 
isoscalar spectral function deduced by analytic continuation from tcN, tctt and NN ^ vrvr 
amplitudes [^]. Note that there is no such thing as a "sigma meson" in this spectrum 
which is completely determined by the (chiral) dynamics of the interacting vrvr and ttN 
system. The integral over ris{t)t~'^ is proportional to the mean squared scalar radius of 
the nucleon. One finds 

(rs)'/'^1.3/m, (10) 



the largest of all nucleon radii, considerably larger than the proton charge radius of 0.86 
fm. 
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2.3. Scalar two-pion exchange 

By its magnitude and range, the form factor Gs{q^) imphes that the nucleon, sur- 
rounded with its two-pion cloud, is the source of a strong scalar-isoscalar field with a 
large effective coupling constant gs = Gs{q^ = 0) = a^jraq ~ 10. When a second 
nucleon couples to this scalar field (see Fig. lb), the resulting two-pion exchange A^A^ 
interaction Vi-k is reminiscent of a Van der Waals force. More than half of the strength 
of V^-K is actually governed by the large spin-isospin polarisability of the nucleon related 
to the transition iV — A in the intermediate state. At long and intermediate distances it 
behaves as 

^2.(r) ~ ^-^P(m^r) , (11) 

where P is a polynomial in niT^r. In the chiral limit (m,r — * 0), this V^-^ approaches the 
characteristic dependence of a non-relativistic Van der Waals potential. 

The two-pion exchange force is the major source of intermediate range attraction that 
binds nuclei. This is, of course, not a new observation. For example, the important role 
of the second-order tensor force from iterated pion exchange had been emphasised long 
ago [ HH] , as well as the close connection of the nuclear force to the strong spin-isospin 
polarisability of the nucleon [I^Ol- The new element that has entered the discussion more 
recently is the systematics provided by chiral effective field theory in dealing with these 
phenomena. 



3. CHIRAL DYNAMICS AND NUCLEAR MATTER 

3.1. Scales at work 

In nuclear matter, the relevant momentum scale is the Fermi momentum kp. Around 
the empirical saturation point with k^p ~ 0.26 GeV ~ 2m^, the Fermi momentum and 
the pion mass are scales of comparable magnitude. This implies that at the densities of 
interest in nuclear physics, p ~ Po = 2(A;^'*)^/37r^ ^ 0.16 fm~'^ ~ 0.45 m^, pions must be 
included as explicit degrees of freedom: their propagation in matter is "resolved" at the 
relevant momentum scales around the Fermi momentum. 

At the same time, kp and m-n are small compared to the characteristic chiral scale, 
47r/^ ^1.2 GeV. Consequently, methods of chiral perturbation theory are expected to be 
applicable to nuclear matter at least in a certain window around kp . In that range, the 
energy density 



A 



p. (12) 



should then be given as a convergent power series in the Fermi momentum. This is our 
working hypothesis. More precisely, the energy per particle has an expansion 

^^^^^ ^Y.^n{kplm.)kl. (13) 



A lOMiv 



n>3 



The expansion coefficients Tn are in general non-trivial functions of kpjm.^^ the dimen- 
sionless ratio of the two relevant scales. These functions must obviously not be further 
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expanded. Apart from kp and mjr, a third relevant "small" scale is the mass difference 
SM = M^—M^ ~ 0.3 GeV between the A(1232) and the nucleon. The strong spin-isospin 
transition from the nucleon to the A isobar is therefore to be included as an additional 
important ingredient in nuclear many-body calculations, so that the J-'n become functions 
of both kp/rriT^ and mT^/5M. 

Let us get a first impression of how the separation of scales controls the scattering 
amplitude T for two nucleons (with momenta \f) \ < kp) interacting in the nuclear medium. 
Omitting its detailed operator structure here for brevity, we denote this amplitude as 
r(g2; kp) and assume that it satisfies, at fixed energy, a subtracted dispersion relation 



with subtraction constants T> . Here > is the (spacelike) squared momentum 
transfered in the t-channel. The amplitudes (with index i) are grouped in a hierarchy of 
terms according to the number of Goldstone bosons (pions) exchanged between the two 
nucleons, or any other mechanisms that may occur at short distance (see Fig.5). The 
corresponding spectral functions Q start at a characteristic threshold /x^ in each channel 
(e.g. = 4m^ for two-pion exchange, /i| = 9m^ for three-pion exchange, and so forth). 
The influence of the Pauli principle in intermediate NN states generally leads to a kp- 
dependence of these spectral functions. 

In the ground state of nuclear matter, the external nucleon lines in Fig.5 have momenta 
in the Fermi sea. The exchanged momentum is limited to Q < 2kp, with an average at 
Q ~ kp. Clearly, in all those processes which have spectral functions starting at fif 3> kp, 
the amplitudes in (14) reduce to the approximate form Tj = T^^^ + df Q"^ where dj = 
(I/tt) dtt~'^(i{t) is proportional to the mean square distance over which this interaction 
takes place. These short-distance contributions with dfQ'^<^l are thus represented as 
contact terms (the constants t/°^), corrected by small finite range (derivative) terms. On 
the other hand, processes with < kp obviously require a full treatment of their spectral 
integrals in Eq.(14). For nuclear densities with kp^ ~ 2m^, one- and two-pion exchange 
processes (as well as those involving low-energy particle-hole excitations) must therefore 
be treated explicitly. They govern the long-range interactions at the distance scales d > 1 
fm relevant to the nuclear many-body problem, whereas short-range mechanisms, with 
t-channel spectral functions involving much larger masses, are not resolved in detail at 
nuclear Fermi momentum scales and can be subsumed in contact interactions. This is 
the "separation of scales" argument that makes strategies of chiral effective field theory 
work even for nuclear problems, with the "small" scales {kp,mTr,SM) distinct from the 
"large" ones (471/,^, M^). The filled Fermi sea of nucleons is an important prerequisite for 
applying low-energy perturbative expansions which would not work in the vacuum. 

3.2. In-medium chiral perturbation theory 

The chiral effective Lagrangian (3) generates the basic pion-nucleon coupling terms used 
in leading-order ChPT. It is of interest [ El |221 to explore whether and to what extent 
this leading-order chiral dynamics can already produce binding and saturation of nuclear 
matter. The new ingredient in performing calculations at finite density (as compared to 




(14) 
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Figure 5. NN amplitude in chiral ef- 
fective field theory: (a) one-pion ex- 
change, (b) two-pion exchange (in- 
cluding A isobar intermediate states), 
(c) contact terms representing short- 
distance dynamics. 




Figure 6. Energy density from 
in-medium chiral perturbation the- 
ory at three-loop order. Dashed 
lines show pions. Each (solid) nu- 
cleon line means insertion of the 
in-medium propagator (15). Not 
shown is an additional two-loop dia- 
gram involving NN contact interac- 
tions. 



evaluations of scattering processes in vacuum) is the in-medium nucleon propagator. For 
a relativistic nucleon with four-momentum = {po,p) it reads 

i^+M^) [^:ri^, - 2^^(p' - M'^)e{po)eikF - ip |)| . (15) 

The second term is the medium insertion which accounts for the fact that the ground 
state of the system has changed from an "empty" vacuum to a filled Fermi sea of nucleons. 
Diagrams can then be organized systematically in the number of medium insertions, and 
an expansion is performed in leading inverse powers of the nucleon mass, consistently 
with the fciT-expansion. 

Our "inward-bound" strategy [ |22j is now as follows. One starts at large distances^ 
(small kp) and systematically generates the pion-induced correlations between nucleons 
as they develop with decreasing distance (increasing kp). Calculations along these lines [ 
122] have been performed to 3-loop order in the energy density (see Fig. 6), including terms 
up to order kp. These calculations incorporate one- and two-pion exchange processes. 
Hartree contributions from one-pion exchange vanish for a spin-saturated system, and 
the one-pion exchange Fock term (upper left of Fig.6) is small, so the leading effect comes 
from the exchange of two pions, with the second order tensor force (in the upper middle 
diagram of Fig.6) providing the dominant attraction. 

At order kp, the momentum space loop integral with iterated pion exchange encounters 
a divergence which needs to be regularised. This can be done in two equivalent ways. 
Either one introduces a high-momentum cutoff A; or one removes the divergent parts 



^Note however that this scheme cannot be pursued down to extremely low densities where new non- 
perturbative phenomena begin to take over, such as formation of nucleon clusters. 
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using dimensional regularisation and replaces them by a counter term. Both procedures 
are equivalent to introducing a subtraction constant in a corresponding dispersion relation 
approach. In either case one ends up with an adjustable NN contact interaction which 
encodes dynamics at short distances not resolved explicitly in the effective low-energy 
theory. This single contact term is the only free parameter at this stage. Adjusting it 
to the binding energy of nuclear matter, the outcome is quite remarkable [122]: despite 
its simplicity, this schematic one-parameter approach already produces nuclear binding 
and saturation (see Fig. 7). One finds E{pq)/A = —15.3 MeV at an equilibrium density 
Po = 0.178 fm'^ together with a (predicted) compressibility K = 255 MeV (empirical: 
K = 220 ±50 MeV [123)- The extension to asymmetric nuclear matter gives a calculated 
asymmetry energy per nucleon, A = 33.8 MeV with no additional parameter (empirical: 
^ ~ 33 MeV [1^). This result suggests that isospin dependent interactions are to large 
extent already accounted for by chiral pion dynamics, without necessity of introducing 
extra ingredients such as the p meson. 

In the scenario just described, entirely based on just leading order chiral pion-nucleon 
dynamics, the binding and saturation of nuclear matter comes from the balance between 
two mechanisms: (a) second order pion exchange with high-momentum cutoff which gives 
driving attraction proportional to kp ^ p in the energy per particle; (b) the Pauli exclusion 
principle acting on intermediate nucleon states of the two-pion pion exchange processes, 
which generates a repulsive contribution of order kp. This 0{kp) term is finite and 
uniquely determined in a model independent way, free of adjustable parameters, involving 
only known low-energy constants {qa, fn) and masses (m,r5 Mn). This statement remains 
unchanged when higher orders in the small-scale expansion are included. 

It is instructive to examine the expansion (13) in the chiral limit (m^r 0) where the 
basic saturation mechanism is already apparent In this limit the coefficients of the 

kp and kp terms of E/A become: 



where small corrections from the one-pion exchange Fock term and from the expansion of 
the kinetic energy have been dropped. One notes that the model- independent JF4 = 0.115 
fm~^ is remarkably close to the empirical JF4 ~ 0.11 fm~^. A high- momentum cutoff scale 
A in the 0.5 GeV range matches the equivalent contact interaction proportional to p for 
which the empirically required strength is JF3 ~ —0.234 fm~^. 

While these exploratory first steps look quite promising, reproducing bulk properties of 
nuclear (and neutron) matter is just a necessary but by no means sufficient condition to 
deal with the nuclear many-body problem. Chiral dynamics with just pions and nucleons 
gives a surprisingly good start for first orientation, but at that (yet incomplete) level of 
the discussion, details such as the momentum dependence of the nucleon single particle 
potential and the density of states at the Fermi surface do not yet come out satisfactorily 
[ |2ni- A further indicator for still missing degrees of freedom is the fact that chiral 
dynamics restricted to pions and nucleons alone gives the critical temperature for the 
nuclear liquid-gas phase transition [E^ too large by about 20-25%. 

The next important step is to incorporate the strong ^ A transition in the two- 
pion exchange dynamics, together with related three-nucleon forces and Pauli corrections. 




(16) 
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Figure 7. Energy per particle of sym- 
metric nuclear matter as determined by 
a one-parameter in-medium ChPT ap- 
proach [123 (solid curve). The dashed 
curve shows for comparison an example 
of a sophisticated many-body calcula- 
tion [I 



Figure 8. Equation of state of symmetric nu- 
clear matter, showing the isotherms of the pres- 
sure P(p, T) as function of baryon density. The 
calculation [ |2H1 includes pions, nucleons and 
A isobars at three-loop order in the free-energy 
density, plus corrections from three-body terms 
at 0{kp) in the energy per particle. 



Such effects were also found to be important in the ab-initio calculations of the Illinois 
group [I2ni- This extended chiral approach [i28i has two subtraction constants encoding 
unresolved short- distance dynamics in the functions J-'s^^ikF/rriTr] m-,^/5M) and a genuine 
three-body term which corrects JFg, but no changes appear in the JF4 term from the 
inclusion of A isobars. The results for both nuclear and neutron matter are now much 
improved. In particular, the momentum dependent (complex) nucleon single-particle 
potential comes out very well in the momentum range p < 400 MeV/c. The Hugenholz - 
Van Hove theorem is strictly fulfilled throughout the procedure. 

In this scenario, nuclear many-body dynamics emerges primarily from the Van der 
Waals - like nature of the attractive two-pion exchange forces in combination with the 
repulsive action of the Pauli principle. It is therefore perhaps not surprising that the 
resulting equation of state at finite temperature [ (5^ (see Fig. 8) is reminiscent of the 
one familiar from a Van der Waals gas. The predicted liquid-gas transition temperature 
Tc ~ 15 MeV is now close to the empirically deduced value, Tc = 16.6 ± 0.9 MeV [^0]. 

At that stage the present approach, guided by chiral effective field theory represent- 
ing low-energy QCD, treats nucleons basically as non-relativistic particles: the low- 
momentum, small-scale expansion of the energy density is simultaneously an expansion in 
inverse powers of the nucleon mass. Contact with the phenomenological non-relativistic 
Skyrme-type energy density functional can be made at this point [I2H1EI1- General fea- 
tures of nuclear Skyrme phenomenology - including surface (gradient) terms - are well 
reproduced around the densities relevant for nuclei (although the chiral approach predicts 
pronounced density dependence of the Skyrme parameters at low densities, induced by 
strong pionic effects when the Fermi momentum becomes comparable to the pion mass). 
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A further interesting comparison can be made with the universal low-momentum NN 
interaction, V^(low-k), derived from a variety of phase-shift equivalent potentials using 
renormalisation group methods ['321- At the Hartree level (linear in the density p), the 
single particle potential generated by V{low-k) agrees well with the one derived from 
chiral dynamics, whereas non-trivial density dependence takes over in the chiral potential 
already at relatively low density, reflecting the action of the Pauli principle on in-medium 
two-pion exchange processes [I2H1- 

There is one place, however, where the non-relativistically reduced chiral theory fails 
completely: it misses the strength of the nuclear spin-orbit force by a large factor. The 
detailed investigation of mechanisms which generate a spin-orbit interaction from two- 
pion exchange (both from the second order tensor force and from terms with intermediate 
A's) reveals that these contributions cancel to a large extent [I2HI- This is a hint that 
genuine relativistic effects, not visible in bulk properties of infinite homogeneous nuclear 
matter, must be considered. 

3.3. Mean fields from QCD condensates 

In relativistic nuclear models [ ^1 1^ , strong Lorentz scalar and vector mean fields, 
each several hundred MeV in magnitude, are at the origin of the abnormally large spin- 
orbit splitting observed in nuclei. While these scalar and vector fields roughly cancel in 
the average single particle potential so that their individual strengths are not revealed 
in the energy per nucleon, they act coherently in building up the spin-orbit potential. 
In-medium QCD sum rules [ give some guidance as to where such strong fields 

have their sources. 

Consider the scalar and vector self-energies Hsyip) appearing in the Dirac equation 
[y^^id^ — Sy) — (Mjv + S5)]\I'Ar = of a nucleon which propagates in isospin-symmetric 
nuclear matter. QCD sum rules establish a relationship between the leading terms of 
and Sv = Sy~°, and the changes with density of the lowest- dimensional quark conden- 
sates, (qq) and {q'^oq) = {q^q) = 'ipj'^- In leading order the condensate part of S5 is 
expressed in terms of the density dependent chiral condensate as follows [ES1E1|: 

= -^lte>, - = ^^^Ps = -^P. . (17) 

with the nucleon scalar density ps = {'^n'^n)- The difference between the vacuum 
condensate {qq)Q and the one at finite density involves the nucleon sigma term (5). The 
Borel scale A^ which roughly separates perturbative and non-perturbative domains in the 
QCD sum rule analysis, has been eliminated in the last step, using the loffe relation (7) 
together with {rriu + rnd){qq) = —m'^f^. The interpretation of Eq.(17) is as follows: once 
nucleons are present in the QCD vacuum, they polarise this vacuum in such a way that 
the magnitude of the (negative) chiral condensate gets reduced with increasing density. 
The resulting attractive scalar mean field is large: 

E5 = MUp) - ~ -6.9 a^— , (18) 

Po 

or S5 ^ —350 MeV at po = 0.16 fm~'^ (using cxn ~ 50 MeV), implying that the nucleon 
mass in nuclear matter is effectively reduced by more than 1/3 of its vacuum value. To the 
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same order in the condensates with lowest dimension, the time component of the isoscalar 
vector self-energy is 

= 3A| ^^'^ = 1^' = <g ' ■ ^''^ 

It reflects the repulsive density-density correlations associated with the time component 
of the quark current, q'j^q. Note that, as pointed out in ref.[|^, the ratio 

^ = Pi (20) 

Ey 4(m„ + nid) p 

is approximately equal to —1 for typical values of the nucleon sigma term and the 
current quark masses m^^rf, and around nuclear matter saturation density where ps — p 
(as an example, take ~ 50 MeV and m„ + ~ 12 MeV at a renormalisation 
scale corresponding to the chiral "gap" 47r/^ ~ 1 GeV). The individually strong scalar 
and vector fields thus balance each other in the average single particle potential, U ~ 
E5 + Sy, and the resulting small contribution simply gets absorbed in the subtraction 
constant attached to the kp term of the chiral expansion of the energy per nucleon. 
However, when introducing local densities for finite systems so that Tis{ps{f)) = S{r) 
and Sy(p(r)) = V{r), the spin-orbit potential 

has a huge strength proportional to ^(0) — 5(0) ~ 0.7 GeV, just the one required by 
phenomenology. 

While useful for further orientation, the QCD sum rule constraints (17-20) are of course 
not very precise at a quantitative level. The estimated error in the ratio S^/Ey ~ — 1 is 
about 20%, given the limited accuracy in the values of ajy and m„ + m^. The leading- 
order loffe formula (7) on which Eq.(17) relies has corrections from condensates of higher 
dimension. In previous QCD sum rule studies the largest uncertainty came from the 
unknown density dependence of the four-quark condensates {qqqq). This uncertainty is 
now considerably reduced, however, by the explicit calculation of chiral two-pion exchange 
contributions to the in-medium nucleon self-energy. 



4. FINITE NUCLEI 

The translation to finite nuclei is best performed using the (relativistic) density func- 
tional framework (for recent reviews, see ref.[|HHl)- Schematically, the energy of the nucleus 
is written as 

E[p] = Tkin + J d^x [Sbgip) + SM] + E,oui . (22) 

with the density p = Y.k=i\4'k){4'k\ expressed in terms of the self-consistent nucleon 
wave functions ipk of the occupied orbits. We associate the "background" part Sbg of 
the energy density with the strong scalar and vector mean fields generated by in-medium 
changes of QCD condensates, while the "exchange correlation" part Sex includes the pionic 
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fluctuations and Pauli principle effects, calculated using in-medium cliiral perturbation 
theory as described in section 3.2. At a practical level, the calculations generating the 
energy density functional are conveniently done by introducing an equivalent effective 
Lagrangian with density dependent four-point couplings between nucleons, 

£e// = ^(^7 ■ d - Mm)^ -W Gi{p) {^Ti^y -W D,{p) (S^T,^)' + C^.^. • (23) 

i i 

While this auxiliary Lagrangian is formally designed to be used in the mean-field limit, 
its physics content reaches out much further. Quantum fluctuations beyond mean field 
are incorporated through the density dependent coupling strengths, Gi{p) and Di{p), 
which multiply the interaction terms {i = S, V, ...) with operators Ts = 1, Ty = etc. 
The matching to the nuclear matter calculations is done at the level of the nucleon self- 
energies expanded in powers of the Fermi momentum, or equivalently, in powers of p^^^. 
This procedure involves the so-called rearrangement terms that occur, as a consequence 
of the density dependent couplings, in order to maintain thermodynamic consistency. 
Electromagnetic interactions (Ce.m.) are also included. 
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Figure 9. Single particle levels for neu- 
trons (left) and protons (right) in ^°Ca 
in comparison with levels deduced from 
experiment. The calculations are per- 
formed in the relativistic point coupling 
model constrained by QCD and chiral 
symmetry [IH7]. 



Figure 10. Single particle levels for protons 
(left) and neutrons (right) in ^^Ni. Calcu- 
lations using Eq.(23) with constraints from 
QCD sum rules and chiral pion dynamics 
(FKVW [EZl) are compared with the rela- 
tivistic mean field model (NL3) of ref . [ I39j . 



Detailed calculations along these lines are reported in [EEIEZI, and the framework is 
summarised and discussed in [ EH] • The self-consistent Dirac equation is solved for the 
single-particle orbits of each given nucleus, so far covering a broad range from ^^O to ^''^Pb 
and to be further extended. The optimal reproduction of binding energies, radii and form 
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factors of the whole set of nuclei requires fine-tuning of the density dependent coupling 
strengths Gi{p). It is remarkable, however, that this fine-tuning systematically remains 
within less than 10% of the constraints set by low-energy QCD in terms of chiral pion 
dynamics and density dependent condensates. The additional surface (derivative) term in 
Eq.(23) turns out to be perfectly consistent with the corresponding gradient terms found 
in the in-medium ChPT calculation of inhomogeneous nuclear matter [ I3H I28j. 

We now present a few representative examples of results from such calculations [I37j. 
Fig. 9 shows the single particle spectra for neutrons and protons in ^^Ca. While the 
binding of these levels is governed by chiral two-pion exchange in the presence of the filled 
Fermi sea (the "Van der Waals plus Pauli" mechanisms that we referred to earlier), the 
spin-orbit splitting is completely determined by the strong scalar and vector mean fields 
that can be interpreted as arising from in-medium changes of the condensate structure of 
the QCD vacuum. Fig. 10 shows our calculations [|37j for ^^Ni in comparison with results 




Figure 11. Charge density distributions and form factors of Ca and Ni calculated 
in the relativistic point coupling model (23) of ref.[EZ] (solid lines). The dashed lines 
are the empirical form factors and charge densities deduced from electron scattering 
experiments. 
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from ref. [ISSl chosen as a representative of the currently most advanced relativistic mean 
field phenomenolgy. Examples of charge distributions and form factors are displayed in 
Fig. 11. These calculations include center-of-mass corrections and convolution with the 
empirical proton form factor. The direct comparison with experimental data is evidently 
quite successful, considering that the detailed surface structure of these distributions is 
reproduced using the gradient terms derived from chiral pion dynamics, with no additional 
fine-tuning required. 

It is instructive to summarise the trends in the equivalent scalar and vector self-energies 
for symmetric nuclear matter, reconstructed using the optimised coupling strengths Gj(p) 
found in the overall "best fit" to finite nuclei. These self-energies have the following 
approximate pattern (with po = 0.16 fm~^): 



(P) 




75 MeV 



0.6 




1/3 



-0.2 




2/3 



0.1 



(24) 



where Jlf ~ -0.34 GeV and ~ +0.34 GeV just balance to zero in their summed 
contribution to the single particle potential but act coherently in producing the large 
spin-orbit splitting. These pieces are compatible with the scalar and vector background 
fields derived from in-medium QCD sum rules. The additional part linear in the density 
p results from contact terms (subtraction constants) in the in-medium ChPT calculation. 
The terms proportional to p^^^ and p^^^ are directly deduced from chiral pion dynamics 
without tuning, while the correction of order p^, representing genuine three-body effects 
and higher order contributions, is subject to adjustment. 

Although the expansion (24) is perturbative by construction, it is a non-trivial feature 
that the detailed density dependence of the nucleon scalar and vector self-energies as 
obtained in ref . [ |3Z| follows very closely, over a wide range of densities, the results of 
Dirac-Brueckner G-matrix calculations [EHl which start from a realistic boson-exchange 
NN interaction. The reasoning behind this observation can presumably be traced to the 
separation of scales at work, a key element of chiral effective field theory which motivates 
the present approach. One- and two-pion exchange in-medium dynamics at momentum 
scales around kp are treated explicitly and include all terms (ladders and others) to 
three-loop order in the energy density. The non-trivial kp dependence in the self-energies 
(beyond "trivial" order kp) reflects in large part the action of the Pauli principle in 
these processes. Such features are also present in the Brueckner ladder summation which 
includes, for example, Pauli blocking effects on iterated one-pion exchange. These model- 
independent terms produce the characteristic kp (or equivalently, p^^^) behaviour in the 
energy per particle. On the other hand, the iteration of the short- distance pieces of the 
NN potential to all orders in the Brueckner ladder, involves intermediate momenta much 
larger than kp and gets absorbed in contact terms which produce self-energy pieces linear 
in the density. Note that these contact interactions must not be iterated further since 
they already represent the full short-distance T-matrix information. 



5. SUMMARY AND OUTLOOK 



The aim of this presentation has been to demonstrate the key role in the structure of 
nucleons and nuclei played by a guiding principle that governs low-energy QCD: spon- 
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taneous chiral symmetry breaking. The effective field theory based on this principle ties 
together a wide range of strong-interaction phenomena involving the lightest {u— and 
(i—) quarks: from the low-energy interactions of pions and nucleons via the condensate 
structure of the QCD vacuum to basic aspects of nuclear binding and saturation. 

We are on the way to a relativistic nuclear energy density functional constrained by 
low-energy QCD. At the least, these constraints significantly reduce the freedom in the 
choice of parameters. In this approach, nuclear binding and saturation arises from the Van 
der Waals - nature of chiral two-pion exchange forces combined with the Pauli principle. 
Strong scalar and vector mean fields, expected to emerge from in-medium changes of QCD 
condensates, largely cancel in the energy per nucleon but drive the spin-orbit splitting in 
finite nuclei. 

The genuine isospin dependence of interactions produced by chiral pion-nucleon dy- 
namics offers interesting perspectives for extrapolations to nuclear systems with extreme 
isospins. A further step, presently taken is the extension of this framework to chiral 
SU(3) dynamics and its application to hypernuclei. 
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